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Abstract —Visible light communication (VLC) has become a 
promising research topic in recent years, and finds its wide 
applications in indoor environments. Particularly, for location 
based services (LBS), visible light also provides a practical 
solution for indoor positioning. Multipath-induced dispersion is 
one of the major concerns for complex indoor environments. It 
affects not only the communication performance but also the 
positioning accuracy. In this paper, we investigate the impact 
of multipath reflections on the positioning accuracy of indoor 
VLC positioning systems. Combined Deterministic and Modified 
Monte Carlo (CDMMC) approach is applied to estimate the 
channel impulse response considering multipath reflections. Since 
the received signal strength (RSS) information is used for the 
positioning algorithm, the power distribution from one trans¬ 
mitter in a typical room configuration is first calculated. Then, 
the positioning accuracy in terms of root mean square error is 
obtained and analyzed. 

Index Terms —Visible light communication. Indoor positioning, 
multipath reflections, impulse response, received signal strength. 

I. Introduction 

IGHT-EMITTING-DIODES (LED) technology has been 
developing rapidly in recent years. It provides efficient 
and economical illumination as well as long service life time. 
Meanwhile, LED technology also finds its wide application 
in visible light communication (VLC) area as LEDs can 
be modulated at relatively high speeds IIl-El. Eor indoor 
environment. Global Positioning System (GPS) signal, serving 
for outdoor positioning, suffers from large attenuation when 
penetrates solid walls. VLC provides a practical solution to re¬ 
alize indoor positioning so that location based services (LBS) 
are available. These systems further offer significant technical 
and operational advantages. Lirst, positioning systems based 
on VLC can be applied in many Radio Lrequency (RE) 
sensitive environments as it does not add any electromagnetic 
interference. Second, employing the same infrastructure for 
both illumination and positioning considerably reduces devel¬ 
opment costs for industry. 

Lor current indoor light positioning systems, LED bulbs act 
as the transmitters, and photo-diode (PD) collocated with a 
user is the receiver. Several positioning approaches have been 
proposed in the literature to calculate the receiver coordinates. 
In one approach, angulation algorithm is used to calculate the 
receiver position based on angle-of-arrival (AOA) information 


and rotation matrix Q. Scene analysis is another approach 
to obtain the receiver position. Leatures of each location are 
collected as the fingerprints in the offline stage. In the online 
stage, the features of current location are measured and by 
matching those with the offline fingerprints, the location of the 
receiver is estimated l8l. In addition, Zigbee technology can be 
combined with VLC to realize long distance positioning ||9l . In 
this paper, we employ a commonly used algorithm where the 
received signal strength (RSS) information is first detected by 
PD, and then the distance between the transmitter and receiver 
is calculated. The lateration algorithm is finally applied to 
estimate the receiver coordinates m, HI. Gaussian mixture 
sigma point particle filter can be further employed to achieve 
high positioning accuracy and prevention of large deviations 

113 . 

In the literature, line-of-sight (LOS) channels have been 
considered without taking into account the multipath reflec¬ 
tions. However, the transmitted signal introduces multipath 
reflections as it bounces off the walls, ceiling and floor where 
the transmitter is a wide-beam LED source, and the receiver 
having a finite field-of-view (LOV) captures reflected photons 
from room surfaces. In this paper, we investigate the effect 
of multipath-induced distortion on the positioning accuracy of 
indoor light positioning systems. There are several methods to 
approximate the impulse response of an indoor optical wireless 
channel. In 03, Barry et al. proposed a deterministic algo¬ 
rithm that partitions a room into many elementary reflectors 
and sums up the impulse response contributions from different 
orders of bounces. As this method is recursive, it takes a 
significant amount of computation time. 

Monte Carlo ray tracing approach is an alternative way 
to calculate the impulse response where rays with identical 
optical power are traced from the source llT4ll . The direction 
of rays are generated by a probability density function (PDL) 
modeled by a Lambertian pattern. When the rays hit reflecting 
surfaces, new rays with reduced power are generated from the 
impact point with the same PDL. This method suffers from 
the fact that it needs a very large number of rays while only 
a small portion of rays will finally reach the PD. To alleviate 
this issue, modified Monte Carlo (MMC) approach is proposed 
which exploits each ray several times instead of only once 
ca. Although this method is fast, it introduces some variance 
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TABLE I 

SYSTEM PARAMETERS 



Fig. 1. System Configuration. 


due to the random direction of the rays. In lIThll . Alqudah 
and Kavehrad proposed a new approach to characterize diffuse 
links in a multiple-input multiple-output (MIMO) system. In 
this paper, we use the methodology developed recently in ini 
referred as combined deterministic and modified Monte Carlo 
(CDMMC), taking advantages of both methods to simulate the 
impulse responses of indoor optical wireless channels. 

The rest of the paper is organized as follows. In Section II, 
the system model and CDMMC approach are briefly discussed. 
In Section III, positioning algorithm is described. In Section 
IV, we present computed numerical results of channel impulse 
response, power distribution and positioning accuracy. Finally, 
Section V concludes the paper. 

II. Multipath Analysis Approach 
A. System Model 

Fig. 1 shows the model of a typical 6 m x 6 m x 3.5 
m indoor optical positioning system under consideration. We 
assume four LED bulbs are installed on the ceiling of a 
room, each of them acts as a transmitter. Each bulb has 
an identification (ID) code assigned to it which denotes its 
coordinates. These LED bulbs are modulated by the drivers 
with on-off-keying (OOK) format. As the bulbs are facing 
downwards, the elevation angle is -90° and azimuth angle is 
0°. Considering the installation, the height of the bulbs is 3.3 
m. The Lambertian mode is order 1 and the wavelength is 420 
nm. A receiver is located at the height of 1.2 m, with an area 
of 1X 10“^m^ and field of view (EOV) of 70°. As it is facing 
upwards, the elevation angle is h- 90° and the azimuth angle 
is 0°. The corresponding reflection coefficients of the walls, 
ceiling and floor along with other parameters are summarized 
in Table 1. 


Room dimensions 

Reflection coefflcients 

length: 6 m 

Pwall- 0.66 

width: 6 m 

POeiling ■ 0.35 

height: 3.5 m 

PFloor- 0.60 

Transmitters (Sources) 

Receiver 

Wavelength: 420 nm 

Area (A): lxl0-"‘m^ 

Height (H): 3.3 m 

Height (/i): 1.2 m 

Lambertian mode (m): 1 

Elevation: +90° 

Elevation: -90° 

Azimuth: 0° 

Azimuth: 0° 

Coordinates: (2,2) (2,4) (4,2) (4,4) 
Power for ”1”/ ”0”: 5 W/3 W 

EOV (To): 70° 


B. Impulse Response Analysis 

CDMMC method combines the deterministic and the MMC 
methods to take advantage of both. The contribution of first 
reflections to the total impulse response is calculated by the 
deterministic method, and then the contributions of second 
and higher-order reflections to the total impulse response are 
calculated by the MMC method. 

In this method, room surfaces are first divided into small 
square elements, each of which has a size of IxlQ-^m^. 
These elements and the PD are considered as the receivers. 
The received power is then obtained as 

p . (0) _ TT (n\ p (0) /ix 

where H (0) is the channel DC gain ifTSll . and Psource^^^ is 
the power from the LED bulbs provided in Table 1. With this 
step, the received power of LOS link is calculated. 

Second, the small elements are considered as point sources 
which their source power can be calculated as 


P, 


source 


(1) = p. 


,( 0 ) 


received P surf ace- 


( 2 ) 


In Eq. (|2]i, Psurface is the reflection coefficient of walls, ceiling 
or floor provided in Table 1. 

Third, these small elements together with PD act as re¬ 
ceivers again, and their received power is 

Preceivef^^ = H (0) Psource^^^ ■ (3) 


With this step, the received power of first reflections is 
calculated. 

To calculate the contribution of second reflections, 10 ran¬ 
dom rays from each of the elements are generated. The power 
Psource''^^ of each ray is then expressed as 




source 


(2) — p . .(t)o f /lO 

— received Hsurjace/ -‘-vv, 


(4) 


where the PDE of the rays’ directions follows 

/ (a, /3) = (“)■ 

In Eq. 0, m is the Lambertion order, a is the angle between 
the random ray vector shown in Eig. 2 and z-axis, and /3 is 
the angle between x-axis and the vector’s projection on the 
X — V plane. Note that Eq. (|5ll is independent of /3. 

The X — Y plane of the unit vector is the surface plane, 
and the origin point is the point source of each ray. These 






















Fig. 2. Coordinate system for a random ray. 



Fig. 3. Flow diagram of CDMMC approach 

rays propagate through the room, and the contribution of the 
secondary reflections is obtained. When the rays hit the room 
surfaces, the impact points are considered as the new sources. 
New rays of each source are generated with reduced power 
where the PDF of their directions follows Eq. (|2i. In this re¬ 
cursive approach, the contribution from subsequent reflections 
is calculated. A flow diagram of CDMMC algorithm is shown 
in Fig. 3. By summing up the contributions from different 
orders of reflections, the total impulse response of the channel 
is computed. 

III. POSITIONING ALGORITHM 

After calculating impulse response with CDMMC method, 
the received signal strength from transmitter i is obtained as 
Pr^\ i = 1,2,3 and 4. The distance di between transmitter 
and receiver is then estimated according to the following 
equation: 

Pr"'' = IT^^COS™ {cj)) Ts (ip) g (ip) cos (ip) Pt. (6) 
Zirdi 

In particular, Eq. (|6]l shows power attenuation from transmitter 
to receiver with respect to their distance where Pt is the 


transmitted power. In (|6l). Parameter (p is irradiance angle of 
transmitter, which numerically equals to the incident angle ip 
of the receiver and is obtained using cos {'tp) = cos {(p) = 
{H — h) /di, where H — h is the vertical distance between 
transmitter and receiver calculated based on the values listed 
in Table 1. Ts{ip) = 1 is the transmittance function of the 
optical filter in the proposed system, and g^tp) is the gain for 
a compound parabolic concentrator (CPC) expressed as 

0 < '0 < ^'c 

V' > 


<7W= 


(7) 


where n denotes the refractive index assumed as n = 1.5, 
and T'c is the FOV provided in Table 1. Therefore, di can be 
calculated as 


di = 


lAT4Pj)g{^)Pt{H-hY 


TT • P, 


( 8 ) 


Given di, the horizontal distance is estimated as 


Ti = \ d, 
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lAT,PiP)gPiP)Pt{H-hY 
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The receiver coordinates (i.e., {x,y)), can be then obtained 
using following equations: 


' {x - xiY + {y - yiY =rl 
(x - X2Y + {y - y2Y =^2 
{x - xsY + {y- ysf = rl 
Xx-XiY + {y-yY^ =rl 


( 10 ) 


where {xi,yi) are the horizontal coordinates of the i* trans¬ 
mitter decoded from the LED ID. Subtracting the last three 
equations from the first one, the following group of equations 
is obtained: 


( 11 ) 


Eq. (fTTl l can be written in matrix format as AX = B where 
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(14) 


X = [a; y] . 

Applying linear least square estimation approach llT9ll . the final 
coordinates are estimated as 

X = (A^A)"^A^B 


(15) 
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Fig. 4. Impulse response of each order of reflections at location A. 


Fig. 5. Impulse response of each order of reflections at location B. 


IV. SIMULATION AND ANALYSIS 

In order to analyze multipath reflections in a typical room, 
three positions are selected. A (0.01, 0.01) represents the 
corner point of the room where strong reflections and scat¬ 
terings happen; B (3, 0.01) represents the edge point where 
the reflections become medium; C (3, 3) represents the point 
at the center of the room where the effect of reflections is 
weak. 

A. Impulse response analysis 

Fig. 4 shows the impulse response at Location A. The total 
impulse response is composed of the LOS, primary reflections 
and secondary reflections. Impulse response amplitude of the 
first reflections is comparable to that of the LOS inducing 
large positioning errors. Compared with the first reflections, 
the impulse response of secondary reflections has a smaller 
amplitude while contributing more delay due to an elongated 
tail. 

Fig. 5 shows the impulse response at Location B. Impulse 
response amplitude from the first reflections is much smaller 
than that of LOS. For the second reflections, the impulse 
response amplitude is much smaller than the first reflections. 
Thus, at this location, as we will show later, the effect of 
reflections on the positioning accuracy becomes less than 
Location A. 

Fig. 6 shows the impulse response at Location C which is 
the center of the room. Compared to the LOS, the amplitudes 
from both reflections are small and can be ignored. Therefore, 
the effect of reflections on the positioning performance is 
expected to be weak in the central area. 

B. Power intensity distribution analysis 

As RSS information is applied in the positioning algorithm 
to estimate the distance between transmitter and receiver, the 
power intensity distribution affects the positioning accuracy 
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Fig. 6. Impulse response of each order of reflections at location C. 


directly. Time division multi-access (TDM) is assumed as 
channel access method; therefore, the four LED bulbs transmit 
signals at different time slots in one frame period Il20ll . In this 
way, PD detects the power intensity from only one transmitter 
at every moment. Considering the symmetrical installation of 
four LED transmitters (TXs), power intensity distribution of 
the bulb located at (2, 2) is estimated where the transmitted 
power is 5 W. The received power of each order is shown in 
Eig. 7 through Fig. 9. 

Particularly, Fig. 7 shows the power intensity distribution 
from LOS. The location which is right below the LED 
bulb has the highest received power. The received power at 
each location is inversely proportional to the distance-squared 
from the source which is in correspondence with the power 
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Fig. 7. Intensity distribution of LOS. 


Fig. 9. Intensity distribution of second reflections. 



Fig. 8. Intensity distribution of first reflections. 


Fig. 10. Positioning error without considering the reflections. 


attenuation expressed as Eq. ©. 

Fig. 8 shows the power intensity distribution contributed 
by the first reflections. The power intensity is high at the edge 
and corner area near the transmitter’s side where the reflections 
are strong. The reflected power will induce positioning errors 
since in the distance estimation, i.e., Eq. (HJ, only direct power 
attenuation from the transmitter is considered. 

Fig. 9 shows the power intensity distribution contributed 
by the second reflections. In the edge and corner area, the 
intensity is slightly higher while the total power intensity 
distribution becomes more uniform. Although the received 
power decreases a lot, these reflections still create positioning 
errors. 


C. Analysis of positioning accuracy 

With the RSS information and positioning algorithm, the 
receiver at each location is able to estimate its coordinates. 
As a benchmark and in order to show the effect of multipath 
reflections on the positioning accuracy, positioning error ne¬ 
glecting the reflected power is also calculated and shown in 
Fig. 10. As it can be seen, the positioning error is low all over 
the room, and only a little higher in the corner area because 
of the low signal-to-noise ratio (SNR). 

Fig. 11, on the other hand, shows the positioning perfor¬ 
mance considering the multipath reflections. The positioning 
performance is satisfactory within the rectangular shown in 
Fig. 11 where the LED bulbs are located right above the 
corners. For the central area, the positioning error is low, due 
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Fig. 11. Positioning error considering the refiections. 


to the small amplitude of impulse response as well as low 
power intensity contributed from reflections. The edge area has 
relatively higher error as the effect of the multipath reflections 
increases. The system has the worst performance at the corner 
area as the reflections are strongest there. 

Fig. 12 and Fig. 13 present the histogram of the positioning 
error of the systems neglecting and considering the reflections, 
respectively. If no reflections are considered, the errors only 
come from the thermal noise and shot noise im which 
are small. The reflections are major concerns in the system 
impairing significantly the system performance. 

Finally, Table 2 compares the positioning error for the two 
cases of neglecting and considering the reflections. When 
reflections are neglected, the root mean square (RMS) error 
is 0.0423 m, and 0.1037 in the corner point. However, this 
situation is not practical as reflections always exist on the 
walls, ceiling and floor in a typical room. 

On the other hand, when multipath reflections are taken 
into account, the positioning accuracy decreases differently 
in distinct locations. At Location A and B, the positioning 
errors are dramatically increased to 2.0989 m and 1.2772 m, 
respectively. As expected, the positioning error at location B 
is less than location A since the reflections at the edge are not 
as strong as those at the corner. However, at the central point 
(i.e.. Location C), the positioning performance is almost not 
affected by the reflections. For the rectangular area covered 
perfectly by the four LED bulbs, the RMS error is 0.0676 m 
which is satisfactory for indoor custom application. The total 
RMS error is 0.8064 m, as the rectangular area covered by 
the LED bulbs is only 11.1% of the total area. Layout design 
of the LED bulbs should be addressed in the future work to 
increase the average positioning accuracy. 
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Fig. 12. Histogram of positioning error with no reflections. 
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Fig. 13. Histogram of positioning error considering the reflections. 


TABLE II 

Positioning error with/without reflections 



Reflections 
neglected (m) 

Reflections 
included (m) 

Loc. A 

0.1037 

2.0989 

Loc. B 

0.0075 

1.2772 

Loc. C 

0.0017 

0.0019 

RMS error of 
the rectangular area 

0.0031 

0.0676 

RMS error of 
the whole room 

0.0423 

0.8064 


















V. CONCLUSIONS 

In this paper, an indoor visible light positioning system 
taking into account multipath reflections has been investigated 
for a typical room where the impulse response is obtained em¬ 
ploying CDMMC approach. In particular, the impulse response 
of each order has been calculated for three locations represent¬ 
ing the corner, edge and center points. As RSS information 
is applied to calculate the distance between a transmitter 
and receiver, the power intensity distribution contributed by 
each reflection has also been computed and shown. Finally, 
positioning errors in all the locations of the room have been 
estimated and compared with those of previous studies where 
reflection is neglected. We have shown that the corner area 
suffers severely from the reflections, while the impact of 
reflections decrease in the edge areas. It has been also shown 
that the central area is almost unaffected by the multipath 
reflections. For future work, calibration methods should be 
considered to increase the positioning accuracy in the corner 
and edge areas, taking into account the multipath reflections. 
In addition, by increasing the coverage area of the LEDs with 
a suitable layout design, the average positioning error can be 
decreased. 
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